Vectorized Solution of ODEs in Matlab

نویسنده

  • Lawrence F. Shampine
چکیده

Vectorization is very important to the efficiency of computation in the popular problem-solving environment Matlab. It is shown that a class of Runge-Kutta methods investigated by Milne and Rosser that compute a block of new values at each step are well-suited to vectorization. Local error estimates and continuous extensions that require no additional function evaluations are derived. A (7,8) pair is derived and implemented in a program BV78 that is shown to perform quite well when compared to the well-known Matlab ODE solver ode45 which is based on a (4,5) pair.

برای دانلود رایگان متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

Matlab : Numerical Solution of Ordinary Differential Equations

Matlab has facilities for the numerical solution of ordinary differential equations (ODEs) of any order. In this document we first consider the solution of a first order ODE. Higher order ODEs can be solved using the same methods, with the higher order equations first having to be reformulated as a system of first order equations. Techniques for solving the first order and second order equation...

متن کامل

Numerical Solution of Delay Differential Equations

After some introductory examples, this chapter considers some of the ways that delay differential equations (DDEs) differ from ordinary differential equations (ODEs). It then discusses numerical methods for DDEs and in particular, how the Runge–Kutta methods that are so popular for ODEs can be extended to DDEs. The treatment of these topics is complete, but it is necessarily brief, so it would ...

متن کامل

An efficient method for the numerical solution of Helmholtz type general two point boundary value problems in ODEs

In this article, we propose and analyze a computational method for numerical solution of general two point boundary value problems. Method is tested on problems to ensure the computational eciency. We have compared numerical results with results obtained by other method in literature. We conclude that propose method is computationally ecient and eective.

متن کامل

picoGA: A Vectorized Implementation of a Genetic Algorithm in Matlab

This document describes picoGA, our submission to the Tiny GA programming competition at GECCO 2006. picoGA is optimized to have small code size and a small memory footprint. It is implemented in the Matlab programming language. Matlab is an interpretted language that is optimized for performing operations on arrays. Loops, especially nested loops, tend to run slowly in Matlab. It is possible h...

متن کامل

Singular boundary value problems for ODEs

This paper is concerned with the numerical solution of a system of ordinary di erential equations ODEs y Sy t f t y p on an inter val b subject to boundary conditions g y y b p The ODEs have a coe cient that is singular at t but it is assumed that the boundary value problem BVP has a smooth solution Some popular methods for BVPs evaluate the ODEs at t This paper deals with the practical issues ...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

عنوان ژورنال:
  • Scalable Computing: Practice and Experience

دوره 10  شماره 

صفحات  -

تاریخ انتشار 2009